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A model is proposed of a collapsing quasi-spherical radiating star with matter content as shear- 
free isotropic fluid undergoing radial heat-flow with outgoing radiation. To describe the radiation 
of the system, we have considered both plane symmetric and spherical Vaidya solutions. Physical 
conditions and thermodynamical relations are studied using local conservation of momentum and 
surface red-shift. We have found that for existence of radiation on the boundary, pressure on the 
boundary is not necessary. 
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I. INTRODUCTION 



In Einstein gravity, gravitational collapse with realistic astronomical matter distribution 
is an important problem of astrophysics. Usually, the formation of compact stellar objects 
such as white dwarf and neutron star are preceded by a period of radiative collapse. Hence 
for astrophysical collapse, it is necessary to describe the appropriate geometry of interior 
and exterior regions and to determine proper junction conditions which allow the matching 
<^ ■ of these regions. 

O 



The study of gravitational collapse was started long ago in 1939 by Oppenheimer 
and Snyder [1]. They studied dust collapse with a static Schwarzschild exterior while 
interior space-time is described by Friedman like solution. Since then several authors have 
extended the above study of collapse of which important and realistic generalizations are 
the following: (i) the static exterior was studied by Misner and Sharp [2] for a perfect fluid 



1 ■ in the interior, (ii) using the idea of outgoing radiation of the collapsing body by Vaidya 



[3] , Santos and collaborations [4-9] included the dissipation in the source by allowing radial 
heat flow (while the body undergoes radiating collapse). Recently, Ghosh and Deskar [10] 
have considered collapse of a radiating star with a plane symmetric boundary (which has 
a close resemblance with spherical symmetry [11]) and have concluded with some general 



remarks. 



So far most of the studies have considered radiating star with interior geometry as spher- 
ical. But in the real astrophysical situation the geometry of the interior of a star may not 
be exactly spherical, rather quasi-spherical in form. So it will be interesting to study quasi- 
spherical interior geometry of a radiating star. In this paper, we have considered the interior 
space-time V~ by Szekeres' model [12, 13] while for exterior geometry V + we have consid- 
ered plane symmetric Vaidya space-time. The plan of the paper is as follows: Exact heat 
flux solution in Szekeres' model has been presented in section II. The junction conditions, 
physical properties and thermodynamical relations are shown in section III. The paper ends 
with a discussion in section IV. 
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II. SOLUTION IN SZEKERES MODEL WITH HEAT FLUX 

The space-time metric for szekeres' model in (n + 2) dimension is in the form [13, 14] 

n 

dsl = -dt 2 + e 2a dr 2 + e 2ti ]T dxj (1) 



i=l 



where a and j3 are functions of all the (n + 2) space-time variables. The stress-energy 
tensor of a non- viscous heat conducting fluid has the expression 

Tpv = (p + p)u^u v + pg^ + q^u v + q v u^ (2) 

where p, p, are the fluid density, isotropic pressure and heat flow vector. We take the 
heat flow vector g M to be orthogonal to the velocity vector i.e., q^u^ = 0. For comoving 
co-ordinate system we choose u M = (1, 0, 0, 0, 0) and q^ = (0, q, qi, Q2, q n ) where q = 
q(t, r, x\, ...,x n ) and qi = qi(t, r, x\, ...x n ), i = 1, 2, n. Now the non- vanishing components 
of the Einstein field equation 

for the above space-time model (1) with matter field in the form of (2) are [13, 14] 

1 ™ f 1 

na{3 + _„(„ - 1)(3 2 - £ j a% + -(n - l)(n - 2)/% + (n - 2)a Xi (3 Xi + a XiXi 



+(« - !)&,*,} + e- 2Q |na'/3' - ^n(n + l)/?' 2 - n/3"| = A + p (3) 

±n(n + l)/3 2 + n/3 - l n (n - l) e - 2 «/?' 2 - e" 2 ^ j±(n - l)(n - 2)/% + („ - l)/3 XiXi j 

= A-p (4) 



d 2 + a + (n - l)d/3 + ^n(n - l)/3 2 + (n - 1)/? + e~ 2 " |(n - l)a'/3' - ^n(n - l)/3' 2 - 
(n - - e- 2 ' 3 V (< + J(« - 2 )( n ~ 3 )^- + + (» " 2 )#w + (« - 3)a x ,/J x A 

- e" 2/3 |(n - l)a Xi /3 Xi + i(„ - l)(„ - 2)/% j = A - p (5) 

(-"xi + Pxi) + Px 3 {ct Xi + (n- 2)0 Xi ) - a XiXj - (n - 2)/3 XiX;J . =0, (i ^ j) (6) 

(/3-d)/3' + /j' = I ge 2« ( 7 ) 



(d - /3)a Xi + d Xi + (n - 1)4, = % e 2/3 (8) 



=0 (9) 
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where dot, dash and subscript stands for partial derivatives with respect to t, r and the 

dp 

<).r, 



corresponding variables respectively (e.g. f3 Xi = ■§§:) and i,j = 1,2,3, ...,n. 



From equations (7) and (9) after differentiating with respect to Xi and t respectively, we 
have the integrability condition 

^ +qa Xi =n[3'p Xi e- 2a , (i = 1, 2, n) (10) 

This equation can not be solved in general. So we have assumed /?' ^ 0, $ Xi = 0. Then 
the form of (3 is 

e = R(t,r)e v( - r ' XuXa -- Xn) (11) 
and from equation (9) we have the solution for a as 

= (12) 
D(t,r) V ; 

where R and D are functions of t, r only. From equations (4) and (5) using equations 
(11) and (12) we have the differential equations of R and D: 

2RR + (n - 1){R 2 — D 2 ) — — (A — p)R 2 = (n - l)f(r) (13) 

n 

and 

RD = f(r) e- 2a (14) 
where f(r) is the arbitrary function of r. 

The function v satisfies the equation 



e^ = A(r)J2^+J2 B ^+ C ^ ( 15 ) 

»=i »=i 

where A, Bi, C are arbitrary functions of r alone with the restriction 



Y,B 2 -AAC = f{r) (16) 
»=i 

Now from equations (7), (8) and (10) we have the components of heat flux vector as 

q=^De~ a (17) 

and 

% = -^^.e-' 3 (18) 

From the above solution, we can see that the field equation (6) is automatically satisfied. 
Using (3), (4) and (5) we have the expression for density as 

•2 , ^ , («+l) , 2A 



(a + n/3 + a 2 + n(3 2 ) + ) -f p + (19) 



(n-l) v " " ' (n-1) * (n-1) 
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However, from equation (14) we note that as R and D are functions of t and r only, so a 
is independent of the space co-ordinates xi 's (i = 1, 2, ...,n) i.e., a Xi = 0, V i = 1,2, ...,n. 
Hence from equation (18) we have gj = and from equation (17) we have seen that 
q = q(t,r) i.e., g is a function of t and r only. Thus only radial heat flow is possible for the 
choice of the metric as we consider. 



III. JUNCTION CONDITIONS AND CONSEQUENCE FOR SZEKERES 
MODEL WITH PLANE SYMMETRIC VAIDYA METRIC 

Let us consider a time-like (n + 1)D hypersurface E, which divides (n + 2)D space-time 
into two distinct (n + 2)D manifolds V~ and V + . For junction conditions we follow the 
modified version of Israel [15] by Santos [4, 5]. Now the geometry of the space-time V~ is 
given by equation (1) while V + and the boundary £ are cheracterised by the metric ansatzs 
as 



~ = 2m(v) 2 _ 2 y 2 

+ (n-l)z"- 1 f-? 



and 



n 



dsl = -dr 2 + A 2 (T)J2<ix 2 (21) 



i=l 



where the arbitrary function m(v) in the Vaidya metric represents the mass at retarded 
time v inside the boundary surface S. Now Israel's junction conditions (as described by 
Santos) are 

(i) The continuity of the line element i.e., 

(ds 2 _)x = (ds 2 + h = dsl (22) 
where ( )s means the value of ( ) on S. 

(ii) The continuity of extrinsic curvature over S gives 

- K+ -K7.= , (23) 
where due to Eiscnhart the extrinsic curvature has the expresion 



jrt = -n±*^-n±r» (24) 



Here ^ l = (r, x±, X2, x n ) are the intrinsic co-ordinates to S, x±i G = 0, 1,2, ...,n+ 1 
are the co-ordinates in V ± and are the components of the normal vector to E in the 
co-ordinates x±- It is to be noted that the above continuity conditions are equivalent to 
junction conditions due to Lichnerowicz and O' Bricn and synge. 

Now for the interior space-time described by the metric (1) the boundary of the interior 
matter distribution (i.e., the surface E) will be characterized by 

/(r,t) = r-r E =0 (25) 

where r s is a constant. As the vector with components is orthogonal to E so we take 



n~ = (0,e a ,0,...,0). 
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So comparing the metric ansatzs given by equations (1) and (21) for dr = we have from 
the continuity relation (22) 



dt 



p$ on T — r (26) 

Also the components of the extrinsic curvature for the interior space-time are 



— = 1, A(t) — e p on r = r v 
dr w 



tf" - and K~ Xi = [P'e 20 -*]^ , i = 1,2, ...,n. 



(27) 



On the other hand for the exterior Vaidya metric described by the equation (20) with its 
exterior boundary, given by 



f(z,v) = 

the unit normal vector to S is given by 



» = 



(28) 



n , 



2m(v) 



-1/2 



(n - l)^™- 1 

and the components of the extrinsic curvature are 



dz 

- — ,1,0,...,0 
dv 



(29) 



d 2 w /du\ 1 m dv 
dr 2 \dr J z n dr 



and 



K+ . = 



dz 
z — 



2m(v) dv 



d,T (n — l)z n 2 dr 



i = 1, 2, n 



(30) 



(31) 



Hence the continuity of the extrinsic curvature due to junction condition (see eq. (23)) 
gives 



and 



2m -(n-2)/3 

(n - 1) dr 



dr J s 



(/3 - /3'e-«) 



(/3 + /3'e- Q ) 



(32) 



(33) 



Now using the junction condition (32) with the help of equation (33), we have 



m(v) = \[{n- l)eS n+1 ^ (/3 2 - /3' 2 er 2a ) 



(34) 



We can intcrprete this as the total energy bounded within the surface S and is equiv- 
alent to the well known mass function in spherical symmetry due to Cahill and McVittie [16]. 

Further using (7), (11), (15), (16), (32), (33) and (34) we obtain (after simplification) 



P s = (qe a h + -n(n-l) [/(r)iT 2 ] E 



(35) 



where the second term on the r.h.s. does not vanish on X. So on the boundary 
vanishing of the isotropic pressure does not imply the vanishing of the heat flux. Thus 



G 



for a quasi-spherical shearing distribution of a collapsing fluid, undergoing dissipa- 
tion in the form of heat flow, the isotropic pressure on the surface of discontinuity S 
does not balance the radiation. Hence in the absence of isotropic pressure there may still 
be radiation on the boundary and the exterior space-time V + will still be Vaidya space-time. 

Moreover, the total luminosity for an observer at rest at infinity is 



Loo = Um (2=i) z n e 
r -> 



/ dm \ 
V dv /£ 



(36) 



If we now consider an observer on the boundary £ then the luminosity for that observer 



is 



is 



n 



dm 
dv 



(37) 



n — 1 



o n0 



(38) 



Thus the boundary red-shift (Z^) of the radiation emitted by a star can be written as 



-PR 



(R + D) 



- 1 



(39) 



Hence the luminosity measured by an observer at rest at infinity is reduced by the 
red-shift in comparison to the luminosity observed on the surface of the collapsing body. 
Also when R + D = then the boundary red-shift attains unlimited value (i.e., Z% — ► oo) 
and the luminosity vanishes at infinity (i.e., L m — > 0). 

We now discuss the thermodynamical relations for a collapsing star. We have seen above 
that when the star produces unpolarizcd radiation while its non-adiabatic fluid collapses 
then the junction condition (35) between pressure p and heat flux q has to be justified. 
Also for physically reasonable fluid we should have (i) < p < p < oo for < r < r s , 
(ii) both & and $Z are negative for r > 0, while < ^ < 1 for < r < r s . Further from 
thermodynamical point of view, we should have the following relations [17]: 

(a) (pzU^).^ = (equation of conservation of matter) where the effective rest mass 
density (measured in the rest frame of u M ), p E is related to the internal energy density U 
by the relation 

p = p E (l + U) 

(b) Gibbs equation: 

TdS = dU + pd(l/p E ) 
where as usual S is the entropy and T is the temperature [18] of the collapsing star. 

(c) Second law of thermodynamics: 



Sft > 
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where the entropy flux 5 M is defined by 

(d) Temperature gradient law: 

<f = -K{g* v + uV)(T, v +Tu v , a u a ) 

with positive thermal conductivity n. But wc note that the second law of thermodynamics 
can be derived from the temperature gradient law, so there is no need to satisfy it. 

Now from the conservation of mass density we have for the present model 

p E =p (r,x u ...,x n ) e -(<"+»« 
with p as effective rest mass density in the infinite past. 

Also from the equation for temperature gradient law the radial heat flow has the form 

BT 

q = -n- e- 2 « (40) 

Further, if we assume the thermal conductivity n as a polynomial in temperature i.e., 

K = 7 T Q > 

then from the comparative study of the expressions for q in equations (17) and (39) we have 
the expression for temperature 

T n+1 = -^±H / (f + ✓) * D)dr + r o(t ) (41) 
with To(t) an arbitrary function of t. 



IV. DISCUSSION 

In this work we have found a general solution for Szekeres' (n + 2)-D space-time model 
with perfect fluid and heat flux. Here we may mention that this solution is not only a 
higher dimensional generalization of Goodc [19] but also it is general in the sense that 
here heat flux is directed along all spatial directions. However, for the study of junction 
conditions with exterior Vaidya model we have considered only radial heat flow. We have 
studied both the physical conditions and the thermodynamical relations for the collapse of 
a radiating star model. 

Moreover, in the Szekeres' model if we consider the co-ordinate transformation [13, 14] 
{xi,x 2 , -,x n ) — > (0i,02, -,0n) by 

x\ = Sin9 n Sin9 n ^i SinQiCot^Qx 

x 2 = Cos6 n Sin6 n _i Sin6 2 Cot^6 1 

xs = Cos8 n -iSind n -2 Sin02Cot^9i 



x n-i = Cos6 3 Sin82Cot^9 1 
x n = Cos02Cot^8i 
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then the form of the (n + 2)-D metric equation (1) becomes 



ds 2 _ = -dt 2 + e 




+ sin 2 9i...sin 2 9 n -idd 2 l ) 
(42) 



Now if we match this interior metric with the spherical (n + 2)-D Vaidya metric i.e., 



which is identical in form to the earlier works of Santos et al [4, 5] and Ghosh et al [11]. 
Also the other conclusions are very much similar to their results. So we have not mentioned 
them here. 

Therefore, for quasi-spherical radiating star with plane symmetric boundary the shear-free 
distribution of collapsing fluid undergoing dissipation in the form of heat flux, it is not 
necessary to have non-vanishing isotropic pressure on the boundary for existence of heat 
flow on it and the result strongly depends on the value of f(r). 
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+sin 2 9i...sin 2 9 n -id9 2 l ) 
(43) 



then on the boundary the relation between pressure and heat flux is 



P E = (9 e ")s 



